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Abstract — This paper considers the so-called multiple-input- 
multiple-output interference channel (MIMO-IC) which has rel- 
evance in applications such as multi-cell coordination in cellular 
networks as well as spectrum sharing in cognitive radio networks 
among others. We consider a beamforming design framework 
based on striking a compromise between beamforming gain at 
the intended receiver (Egoism) and the mitigation of interference 
\ created towards other receivers (Altruism). Combining egoistic 
and altruistic beamforming has been shown previously in several 
' papers to be instrumental to optimizing the rates in a multiple- 
' input-single-output interference channel MISO-IC (i.e. where 
receivers have no interference canceling capability). Here, by 
using the framework of Bayesian games, we shed more light 
on these game-theoretic concepts in the more general context 
of MIMO channels and more particularly when coordinating 
parties only have channel state information (CSI) of channels that 
they can measure directly. This allows us to derive distributed 
beamforming techniques. We draw parallels with existing work 
on the MIMO-IC, including rate-optimizing and interference- 
\ alignment precoding techniques, showing how such techniques 
, may be improved or re-interpreted through a common prism 
• based on balancing egoistic and altruistic beamforming. Our 
analysis and simulations currently limited to single stream trans- 
mission per user attest the improvements over known interference 
alignment based methods in terms of sum rate performance in 
I the case of so-called asymmetric networks. 

Index Terms — multi-cell, MIMO, distributed beamforming, 
Pareto boundary, game theory, Bayesian equilibrium, interfer- 
ence channels, distributed bargaining, egoistic, altruistic, inter- 
, ference alignment 



I. Introduction 

The mitigation of interference in multi-point to multi-point 
radio systems is of utmost importance and has relevance in 
several practical contexts. Among the more popular cases, 
we may cite the optimization of multi cell multiple-input- 
multiple-output (MIMO) systems with full frequency reuse 
and cognitive radio scenarios featuring two or more service 
providers sharing an identical spectrum license on overlapping 
coverage areas. In all these cases, the system may be modeled 
as a network of Nc interfering radio links where each link 
consists of a sender trying to communicate messages to a 
unique receiver in spite of the interference arising from or 
created towards other links. 

For system limitation or privacy reasons, when the backhaul 
network cannot support a complete sharing of data sym- 
bols across all transmitters (Txs), the channel remains an 
interference channel (IC). Coordination in terms of beam- 
forming is required to be decentralized in the sense that 
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global channel state information at transmitters (CSIT) may 
not be available everywhere. In the context of distributed 
beamforming, game theory appears as a sensible approach as 
a basis for algorithm design. Recently an interesting game 
theory framework for beamforming-based coordination was 
proposed for the multiple-input-single-output (MISO) case 
by which the transmitters (e.g. the base stations) seek to 
strike a compromise between selfishly serving their users 
while ignoring the interference effects on the one hand, and 
altruistically minimizing the harm they cause to other non- 
intended receivers on the other hand. An important result in 
this area was the characterization of all so-called Pareto rate 
optimal beamforming solutions for the two-cell case in the 
form of positive linear combinations of the purely selfish and 
purely altruistic beamforming solutions [l]-[3] and [4]-[6] in 
the case of partial CSI. Unfortunately, how or whether at all 
this analysis can be extended to the context of MIMO interfer- 
ence channels (i.e. where receivers have themselves multiple 
antennas and interference cancelling capability) remains an 
open question. 

In parallel, coordination on the MIMO interference channel 
has emerged as a very popular topic in its own right, with 
several important non-game related contributions shedding 
light on rate-scaling optimal precoding strategies based on 
so-called interference alignment, subspace optimization, alter- 
nated maximum signal to interference and noise ratio (SINR) 
optimization, [7]-[9] and rate-maximizing precoding strategies 
[10], [11], to cite just a few examples. 

Interference alignment based strategies exhibit the designed 
feature of rendering interference cancellable (when feasible, 
according to the available degrees of freedom) at both the 
transmitter and receiver side. Such a behaviour is optimal in 
the large signal to noise ratio (SNR) region when Rxs have 
single user decoder and interference is the key bottleneck. At 
finite SNR, various strategies exist which aim at maximizing 
a link quality metric individually over each link, while tak- 
ing interference into account. This often takes the form of 
maximizing the link's SINR or minimizing minimum-mean- 
square-error (MMSE). This approach provides good rates in 
symmetric networks where all links are subject to impairments 
(noise, average interference) of similar level. In more general 
and practical situations however, we argue that a better sum 
rate may be obtained from a proper and different weighting 
of the egoistic and altruistic objective at each individual link. 
This situation is particularly important when more links are 
subject to statistically stronger interference than others, a 
case which has so far received little attention and which we 
shall refer here as asymmetric networks. For this purpose, we 
suggest to re-visit the problem of coordinated beamforming 
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design by directly building on the game theoretic concept of 
egoistic and altruistic game equalibria. Because our focus is 
on scenarios where CSI is not fully available, we consider 
a class of games suitable to the case of partial information- 
based decision making, called Bayesian games. Note that this 
is different from the limited CSI feedback scenario studied 
by previous authors [12] who consider channel quantization 
requirement as function of SNR. Our approach is two fold, first 
derive analytically the game equalibria. Second, exploit the 
obtained equilibria solution into heuristic design of a practical 
beamforming teachnique. The behaviour of our solution is then 
studied both theoretically (large SNR regime) and tested by 
simulations. 

More specifically in this paper, our contributions are as 
follows: 

• We define the egoistic and altruistic objective functions 
and derive analytically the equilibria of so-called egoistic 
and altruistic Bayesian games [13]. 

• Based on the equilibria, we propose a practical distributed 
beamforming scheme which provides a game-theoretic 
interpretation of the distributed sum rate maximization 
problem the MIMO-IC, such as [11]. 

• The proposed techniques allows a tradeoff between the 
reduced complexity/feedback and the rate maximization 
offered by [11]. 

• We show that our algorithm exhibits the same rate 
scaling (when SNR grows) as shown by recent interest- 
ing interference alignment based methods [7]-[9] which 
operates on the same feedback assumption as the pro- 
posed beamforming scheme. At finite SNR, we show 
improvements in terms of sum rate, especially in the case 
of asymmetric networks where interference-alignment 
methods are unable to properly weigh the contributions 
on the different interfering links to maximize the sum 
rate. This situation is particularly relevant. In practical 
contexts where for complexity limitation reasons only 
a subset of cells (links) is coordinated across, while 
other uncoordinated links contribute to additional unequal 
amounts of unstructural interference. 

A. Notations 

The lower case bold face letter represents a vector whereas 
the upper case bold face letter represents a matrix. (.)^ 
represents the complex conjugate transpose. I is the identity 
matrix. T/('"'^=^)(A) (resp. (A)) is the eigenvector cor- 

responding to the largest (resp. smallest) eigenvalue of A. £b 
is the expectation operator over the statistics of the random 
variable B. §\M define a set of elements in § excluding the 
elements in B. Tr{A) denotes the trace of matrix A. 

II. Bayesian Games Definition on interference 

CHANNEL 

Let J\f = {!,..., N} be a set containing a finite set Afc, 
with cardinality Nc < N, of cooperating transmitters (Txs), 
also termed as players. From now on, we use players and 
Txs interchangably. We call the set Afc a coordination cluster 
and Txs outside the cluster will contribute to uncontrolled 



interference. The provided model has general applications in 
which the Txs can be base stations in cellular downlink where 
typically coordination is restricted to a subset of neighbouring 
cell sites while more distant sites cannot be coordinated over 
[14] ; nodes in ad-hoc network and cognitive radio. 

Each Tx is equipped with Nt antennas and the Rx with 
Nr antennas. Each Tx communicates with a unique Rx at 
a time. Txs are not allowed or able to exchange users' 
packet (message) information, giving rise to an interference 
channel over which we seek some form of beamforming-based 
coordination. The channel from Tx i to Rx j Hj,; E QNrxNt 
is given by: 

Hjj = ^/a~Hji, i, j = 1, . . . ,iVc (1) 

Each element in channel matrix Hj^ is an independent 
identically distributed complex Gaussian random variable with 
zero mean and unit variance and aji denotes the slow- 
varying shadowing and pathloss attenuation. Hj^ is circularly 
symmetric complex gaussian and the probability density is 

/h,,(H) = -±^exp{~Tr (HH^)). (2) 

A. Limited Channel knowledge 

Although there may exist various ranges and definitions 
of local CSI, we assume a standard definition of a quasi- 
distributed CSI scenario where the devices (Tx and Rx alike) 
are able to gain knowledge of those local channel coefficients 
directly connected to them, as illustrated in Fig. [T] possibly 
complemented with some Umited non local information (to be 
defined later). 

The set of CSI locally available (resp. not available) at Tx 
i denoted by (resp. B^) is denoted by: 

Bi = ; B^ = {Rki}k,i=i...N^ \ (3) 

Similarly, define the set of channels known (resp. un- 
known) at Rx i denoted by (resp. M^) as: = 

{Hjj}j=i,...,Ar, ; = {Hfei}fc,i=i...Ar, \ ^^■ By Construc- 
tion here, locally available channel knowledge, Bi, is only 
known to Tx i but not other Txs. We call this knowledge Bi 
the type of player (Tx) i, in the game theoretic terminology 
[13]. 

In the view of Tx i, the decision to be made shall be based 
on its type B^ and its beliefs on other Txs types. Since Tx i 
does not know other Txs types, we assume that Tx i has a 
probability density over the possible values of other players 
channel knowledge Bj. For simplicity, we assume that these 
beliefs are symmetric: the probability density of the gaussian 
channels available at Tx i regarding Bj is the same as the 
probability density of Tx j over B^. The asymmetric path loss 
antennuations aji are assumed to be long term satistics and 
known to the Txs. And we assume that the channel coefficients 
in the network are statistically independent from each other 
We define here the joint beliefs (probability density) at Tx i: 

Ai.=p(B^)=/H,,(H)^^^^^-^^=M. (4) 

The Tx index i is dropped because the beliefs are symmetric 
among Txs, given the asymetric path loss coefficients aji. p{.) 
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Txi Tx2 Tx3 Tx4 Txs "^^Nc 



Ml 




Rxi Rx2 Rx^ Rx4 Rxf^ Rx 

Fig. 1: Limited channel knowledge model: as an illustration, 
the local CSI available at Tx Nc is shown in dashed lines. The 
local CSI available at Rx 1 is shown in solid lines. 



is a probability measure and /g^ (H) is the probability density 
of a complex gaussian channel defined in (|2]i. The second 
equality relies on the assumptions that the channel coefficients 
from any Tx to any Rx are independent. 

Based on its belief, Tx i designs the transmit beamforming 
vector, Wj e C^*^^. As in several important contributions 
dealing with coordination on the interference channel [2], [8], 
[15]-[19], we assume linear beamforming. We call the transmit 
beamforming vector Wj an action of Tx i and denote the set 
of all possible actions by A at any Tx. 

yl= {weC^'^i : ||w||2 < 1} (5) 
The received signal at Rx i is therefore 

Di = vf HjiWj + ^ vf + (6) 

where rii is a gaussian noise with power af. Note that the 
noise levels of depend on the link index which was not 
considered in previous work on transmitter coordination. The 
Rxs are assumed to employ maximum SINR (Max-SINR) 
beamforming throughout the paper so as to also maximize 
the link rates [20]. The receive beamformer is classically 
given by: 

Vj = — (7) 

||C_Ri Hj,,Wj|| 

where C^,; is the covariance matrix of received interference 
and noise 

Cffi = Hy w^wf H,^P + afl. (8) 

P is the transmit power Note that the receive beamformer 
is a function of all transmit beamforming vectors w^. When 
the transmit beamforming vector w, is optimized, the received 
beamforming vector is modified accordingly. 




N = 7; N,, = A 

Fig. 2: This figure illustrates a system of = 7 cells where 
A^c = 4 form a coordination cluster Empty squares repre- 
sent transmitters whereas filled squares represent receivers. 
The noise power (which includes out of cluster interference) 
undergone in each cell varies from link to link. 



Importantly, the noise will in practice capture thermal noise 
effects but also any interference originating from the rest of 
the network, i.e. coming from transmitters located beyond 
the coordination cluster. Thus, depending on path loss and 
shadowing effects, the {cr^} may be quite different from each 
other [21]. Fig. |4] illustrates a system of A^ = 7 cells where 
A'c = 4 form a coordination cluster Note that we consider the 
sum of uncoordinated source of interference and thermal noise 
to be spatially white. The non-colored interference assumption 
is justified in the scenario where receivers cannot obtain 
specific knowledge of the interference covariance and can 
be interpreted as a worst case scenario, since the receivers 
cannot use their spatial degrees of freedom to further cancel 
uncontrolled interference. 

Receiver feedback v.s. Reciprocal Channel: In the case 
of reciprocal channels , e.g. time-devision-duplex systems 
(TDD), the feedback requirement to obtain B,; can be replaced 
by a channel estimation step based on uplink pilot sequences. 
Additionally, it will be classically assumed that the receivers 
are able to estimate the covariance matrix of their interference 
signal, based on, say, transmit pilot sequences. 

We can now define the Bayesian game on interference 
channel. 

Definition 1: The Bayesian game on interference channel 
can be described by a 5-tuple: 

G=<AA„A{Bj,M,{wJ > (9) 

where fi denotes the beliefs of the players and {ui] denotes 
the utility functions of the players, which can be either egoistic 
or altruistic. 

Specific definitions of Ui will be given in the following 
sections. The players are assumed to be rational as they 
maximize their own utility based on their types and beliefs. 
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Definition 2: A pure-strategy of player i, Si : Mi is a 

deterministic choice of action given information of player 
i. 

Definition 3: A strategy profile s* — (s*, slj achieves the 
Bayesian Equilibrium if s* is the best response of player i 
given strategy tuple slj for all other players and is character- 
ized by 

Vi s* = argmaxfgi {ui(si,slj)} . (10) 
Note that, intuitively, the player's strategy is optimized by 
averaging over the beliefs (the distribution of all missing state 
information) while in a standard game, such expectation is not 
required. 

In the following sections, we derive the equilibria for 
egoistic and altruistic bayesian games respectively. These 
equilibria constitute extreme strategies which do not perform 
optimally in terms of the overall network performance, yet can 
be exploited as components of a more general beamforming- 
based coordination technique which is then proposed in section 

m 

III. Bayesian Games with Receiver Beamformer 
Feedback 

We assume that Tx has the local channel state information 
Mi and the added knowledge of receive beamformers through a 
feedback channel. Note that in the case of reciprocal channels, 
the receive beamformer feedback is not required. 

A. Egoistic Bayesian Game 

Definition 4: Denote the set of transmit beamforming vec- 
tors of players j,j ^ i, by w„i. The egoistic utility function 
for Tx i is defined as its received SINR 



-Ui(Wi, W_j) 



vfH„,w,;PP 



(11) 



Based on Tx i's belief, Tx i maximizes the utility function in 
( fTTT i where Vj is a known quantity. 

Lemma 1: There exist at least one Bayesian EquiUbrium in 
the egoistic Bayesian Game G Q with utility function defined 
in (El)- 

Proof: Ai is convex, closed and bounded for all players 
i and the egoistic utility function Mi(wi,w_i) is continuous 
in both Wi and w_i. The utility function is convex in for 
any set w_i. Thus, at least one Bayesian Equilibrium exists 
[22], [23]. ■ 
Theorem 1: The best-response strategy of player i in the 
egoistic Bayesian Game G Q with utility function (fTTT i is to 
maximize the utility function based on its beUef: 

wf^° = argmax^B^ {ui(wi, w_i)} . (12) 

The best-response strategy of player i is 

^Ego ^y{,nax)^^^^ (13) 

where denotes the egoistic equilibrium matrix for Tx i, 
given by 



Proof The knowledge of receive beamformers decorre- 
lates the maximization problem which can be written as 



Ego 



arg max £j. 



1 



wfH,^v,vfH,,w, 



(14) 



The egoistic-optimal transmit beamformer is therefore the 
dominant eigenvector of H^Vivf^Hii. ■ 

B. Altruistic Bayesian Game 

Definition 5: The utility of the altruistic game is defined 
here so as to minimize the sum of interference powers caused 
to other receivers. 



(15) 



Note that the receive beamforming vectors Vj is a Max-SINR 
receiver which depends on the transmit beamforming vectors 
•Wj and cause conflicts between Txs. 

Lemma 2: There exist at least one Bayesian Equilibrium 
in the altruistic Bayesian Game G (O with utility function 
defined in ( fTSI ). 

Proof: Ai is convex, closed and bounded for all players 
i and the altruistic utility function Mi(wi,w_,:) is continuous 
in both Wi and w_i. The utility function is concave in w,; for 
any set w_i. Thus, at least one Bayesian Equilibrium exists 
[22], [23] . ■ 

Theorem 2: Based on belief /i, Tx i seeks to maximize the 
utility function defined in ([TST l. The best-response strategy is 



Alt 



(16) 



where A.ji denotes the altruistic equilibrium matrix for Tx i 
towards Rx j, defined by Aji — Hj^v^v^Hji. 

Proof: Recall the utility function to be 



Hj.Wil^ = -Ej5^«wf AjiW,. Since v^- 
are known from feedback or estimation in reciprocal 
channels, the optimal is the least dominant eigenvector of 
the matrix ^ji- ■ 

IV. SuMRATE Maximization with Receive 
Beamformer Feedback 

From the results above, it can be seen that balancing 
altruism and egoism for player i can be done by trading- 
off between setting the beamformer close to the dominant 
eigenvectors of the egoistic equilibrium E, or that of the 
negative altruistic equilibrium {—Aji} (j ^ i) matrices in 
([T6] l. Interestingly, it can be shown that sum rate maximizing 
precoding for the MIMO-IC does exactly that. Thus we hereby 
briefly re-visit rate-maximization approaches such as [11] with 
this perspective. 

Denote the sum rate by ^ = X^iLi where Ri = 

|vf HiiWil^P 



l0g2 1 



5:,';ijvfH.,w,|2p+a? 
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Lemma 3: The transmit beamforming vector which maxi- 
mizes the sum rate R is the dominant eigenvector of a matrix, 
which is a hnear combination of Ei and A,v; 



E,; 



where 



(17) 



(18) 



where S'jfc = jv^Hj^WfepP and /imax is defined in the proof. 

Proof: see appendix IVIII-AI ■ 
Note that the balancing between altruism and egoism in sum 
rate maximization is done using the dominant eigenvector 
of a simple linear combination of the altruistic and egoistic 
equilibrium matrices. The balancing parameters, {A°f*}, can 
be shown simply to coincide with the pricing parameters 
invoked in the iterative algorithm proposed in [11]. Clearly, 
these parameters plays a key role, however their computation is 
a function of the global channel state information and requires 
additional message (price) exchange. Instead, we seek below 
a suboptimal egoism-altruism balancing technique which only 
requires statistical channel information, while exhibiting the 
right performance scaling when SNR grows large. 

V. A PRACTICAL DISTRIBUTED BEAMFORMING 
ALGORITHM: DBA 

We are proposing the following distributed beamforming 
algorithm {DBA) where one computes the transmit and receive 
beamformers iteratively as: 



We 



yn 



Ei + ^ AjjAj 



(19) 



(20) 



where \ji shall be made to depend on channel statistics 
only. At this stage, it is interesting to compare with previous 
schemes based on interference alignment such as the practical 
algorithms proposed in [9]. In such schemes, the transmit 
beamformer is taken independent of 'Ha. Note that here 
however, w, is correlated to the direct channel gain 
through the Egoistic matrix E^ in DBA. The correlation is 
useful in terms of sum rate as it allows proper weighting 
between the contributions of the egoistic and altruistic matrices 
in a link specific manner. 

A. The egoism-altruism balancing parameters Xji 

The egoism-altruism balancing parameters Xji are now 
found heuristically based on the statistical channel informa- 
tion. Recall from ( fTsT i that 



5, 



S, 



(21) 



where S, = |vf H,,w,f P and /, = E.^, |vf H,feW,f P 



Following the principle behind sum rate maximization, we 
conjecture that at convergence, residual coordinated interfer- 
ence shall be proportionate to the noise and out-of-cluster 
interference, i.e. Ij — 0{<j^). Note that this should not be 
interpreted as an assumption in a proof but rather as a proposed 
design guideline. Based on this, we propose the following 
characterization: 



5, 



5. + 0(a2) 



Note that Sj and Si are independent and we have 



A;r* 



exf 



E 



S, + 0(^2) 



> 



£{S.) + Oiaf) 



£S, 



OK) 



0(^1) 



1 



ok: 



(22) 

(23) 
(24) 
(25) 

(26) 



where (a) is because 5,; , Sj are independent and (&) is because 
the function is concave in x and therefore by Jensen's 
inequality, we have J''^. > . 

^ tx+c — x+c 

Although £Si is not known explicitly, it is strongly related 
to the strength of the direct channel Pan. Let 7^ = In 
order to obtain an exploitable formulation for Xji, we replace 
SSi by Pail and 0{af) by af, to derive: 



Xji 



1 



1 + 7- 



1+77 



(27) 



Pan 



Interestingly, in the special case where direct channels have 
the same average strength, we obtain a simple expression 



A, 



l + 7[ 



1 + 7 



— 7j- 



(28) 



The above result suggests Tx i to behave more altruistically 
towards link j when the SNR of link j is high or when the SNR 
of link i is comparatively lower. This is in accordance with 
the intuition behind rate maximization over parallel gaussian 
channels. 

DBA iterates between optimizing the transmit and receive 
beamformers, as summarized in Algorithm [T] Iterating be- 
tween transmit and receive beamformers is reminiscent of 
recent interference-alignment based methods [8], [9]. However 
here, interference alignment is not a design criterion. In [8], 
an improved interference alignment technique based on alter- 
nately maximizing the SINR at both transmitter and receiver 
sides is proposed. In contrast, here the Max-SINR criterion 
is only used at the receiver side. Although the distinction is 
unimportant in the large SNR case (see below), it dramatically 
changes performance in certain situations at finite SNR (see 
Section IVB. 
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Algorithm 1 DBA 



1) Initialize beamforming vectors w^, i = 1, . . . , Nc, to be 
predefined vectors. 

2) For each Rx compute v; — h„w. 



l|C^-H„w,|| 

computed with Wj in previous step. 
3) For each Tx i, compute 

rraax I tt" i v^-^^ 



+ E7#^ ^J'Aj* ) where 



A 



ymax 

from satistical parameters ( |27| i 
4) Repeat step 2 and 3 until convergence 



where Cjn is 
are computed 



B. Asymtotic Interference Alignment 

One important aspect of the algorithm above is whether it 
achieves the interference alignment in high SNR regime [8]. 
The following theorem answers this question positively. 

Definition 6: Define the set of beamforming vectors solu- 
tions in downlink (respectively uplink) interference ahgnment 
to be [8] 



TA 



DL 



{(wi, . . . , watJ 



lA 



UL 



HifeWfeW^H^ is low rank, Vi 
{(vi,. 

E 



■ jVatJ : 



H^VfcV^Hfci is low rank, Vi 



(29) 



(30) 



Thus, for all (w^, , 
beamformers Vi,i 
satisfied: 



, wjv^) e lA^'", there exist receive 
1, . . . ,Nc such that the following is 



vfH,jW, =0 V^,j^^. 



(31) 



Note that the uplink alignment solutions are defined for a 
virtual uplink having the same frequency and only appear here 
as a technical concept helping with the proof. 

Theorem 3: Assume the downlink interference alignment 
set is non-empty (interference alignment is feasible). Denote 
average SNR of link i by 7^ = 
then in the large SNR regime, P 



Let Aj, = -T-^lj 



any transmit 

beamforming vector in 2A^^ is a convergence (stable) point 
of DBA. 

Proof: see Appendix I VIII-BI ■ 
Note that this does not prove global convergence, but local 
convergence, as is the case for other lA or rate maximization 
techniques [8], [9], [11]. Another way to characterize local 
convergence is as follows: assuming interference alignment 
is feasible (2A^^ is non-empty), the first algorithm in [8] 
was shown to converge to transmit beamformers belonging to 
2A^^ and the receivers are based on the minimum eigen- 
vector of the dowlink interference covariance matrix, which 
tends to be low-rank. However, DBA selects its receive beam- 
former from the Max-SINR criterion which, in the large SNR 
situation, is also identical to selecting receive beamformers in 
the null space of the interference covariance matrix. Therefore 
when interference alignment is feasible, the algorithm in [8] 
and DBA coincide at large SNR. This aspect is confirmed by 
our simulations (see section IVII ). 




10 12 

SNR (dB), |N_^,N|,NJ=[3,2,2] 



Fig. 3: Sum rate comparison in multi links systems is il- 
lustrated with [Nc,Nt,Nr] = [3,2,2] with increasing SNR. 
DBA, SR-Max and Max-SINR achieve very close performance 
in symmetric networks. 



VI. Simulation Results 

In this section, we investigate the sum rate performances 
of DBA in comparison with several related methods, namely 
the Max-SINR method [8], the alternated-minimization (Alt- 
Min) method for interference alignment [9] and the sum rate 
optimization method (SR-Max) [11]. The SR-Max method is 
by construction optimal but is more complex and requires 
extra sharing or feedback of pricing information among the 
transmitters. To ensure a fair comparison, all the algorithms 
in comparisons are initialized to the same solution and have 
the same stopping condition. The algorithms are considered to 
reach convergence if the sum rates achieved between succes- 
sive iterations have difference less than 0.001. We perform sum 
rate comparisons in both symmetric channels and asymmetric 
channels where links undergo different levels of out-of-cluster 
noise. Define the Signal to Interference ratio of link i to be 
SIRi ~ ^n" — . The SIB is assumed to be 1 for all links, 

unless otherwise stated. Denote the difference in SNR between 
two links in asymmetric channels by AS NR. Note that the 
proposed algorithm is not limited to the following settings, but 
can be applied to network with arbitrary players and number 
of antennas. 



A. Symmetric Channels 

Fig. [3] illustrates the sum rate comparison of DBA with Max- 
SINR, Alt-Min and SR-Max in a system of 3 links and each 
Tx and Rx have 2 antennas. Since interference alignment is 
feasible in this case, the sum rate performance of SR-Max and 
Max-SINR increases linearly with SNR. DBA achieves sum 
rate performance with the same scaling as Max-SINR and SR- 
Max (i.e. multiplexing gain of 3). Therefore these methods 
seem to perform similarly in symmetric channels. 
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B. Asymmetric Channels 

In the asymmetric system, some links undergo uneven levels 
of noise and uncontrolled interference. Another aspect is that 
more links can experience greater path loss or shadowing than 
others. Here we consider a few typical scenarios for which 
could constitute asymmetric networks, as shown in Fig. |4] 

1) Asymmetric uncontrolled interference power, illustrated 
in Fig. Wa\ In Fig. |5j there are 3 links in the system in which 
the noise and unstructural interference in one of the links are 
20dB stronger than the other two links. This set up captures the 
scenario that one link is at the boundary of the coordination 
cluster and suffer from strong out-of-cluster noise. The SIR 
of every link is assumed to be 10 dB. in this scenario, DBA 
outperforms interference alignment based methods because 
they are unable to properly weigh the importance of each link 
in the overall sum rate. SR-Max is by construction sum rate 
optimal. However, in the asymmetric network, we observe by 
simulation that the convergence may require more iterations 
than other algorithms and the increment in sum rate per 
iteration can be small in some channel realizations. 

2) Asymmetric uncontrolled interference power and inter- 
ference within cluster, illustrated in Fig. \4b\ In Fig. |6l we 
compare the sum rate performance in the same set up as 
in Fig. m except that the SIR's of the Hnks are [10, 10,0.1] 
respectively. Thus, link 3 not only suffers from strong out 
of cluster noise, but also suffers from strong interference 
within the cluster The asymmetry penalizes the Max-SINR 
and interference alignment methods because they are unable 
to properly weigh the contributions of the weaker link in 
the sum rate. The Max-SINR strategy turns out to make link 
3 very egoistic in this example, while its proper behavior 
should be altruistic. In contrast, DBA exploits useful statistical 
information, allowing weaker link to allocate their spatial 
degrees of freedom wisely towards helping stronger links and 
vice versa, yielding a better sum rate for the same feedback 
budget. The performance is very close to SR-Max, with less 
information exchange. 

3) Asymmetric desired channel power, illustrated in Fis. \4c\ 
In Fig. |2l there are 3 links cooperating in the system. Each 
Tx and Rx has 2 antennas and has 1 stream transmission. The 
noise at each Rx is the same. The system is asymmetric in 
a sense that the direct channel gain Hn of link 1 is 30dB 
weaker than other links in the network. This set up models a 
realistic environment where the user suffers strong shadowing. 
DBA achieves sum rate closed to SR-Max and much better 
than other interference alignment based schemes Max-SINR 
and Alt-Min. 

VII. Conclusion 

We model the distributed beamforming optimization prob- 
lem on MIMO interference channel using the framework of 
Bayesian Games which allow players to have imcomplete 
information of the game, in this case the channel state 
information. Based on the incentives of the players, we 
proposed two games: the Egoistic Bayesian Game (players 
selfishly maximize its rate) and the Altruistic Bayesian Game 
(players altruisticly minimize interference generated towards 



other players). We proved the existence of equilibria of such 
games and the best response strategy of players are computed. 
Inspired from the equilibria, a beamforming technique based 
on balancing the egoistic and the altruistic behavior with the 
aim of maximizing the sum rate is proposed. Such beamform- 
ing algorithm exhibits the same optimal rate scaling (when 
SNR grows) shown by recent iterative interference-aUgnment 
based methods. The proposed beamforming algorithm acheives 
close to optimal sum rate maximization method [11] without 
additional pricing feedbacks from users and outperform in- 
terference alignment based methods in terms of sum rate in 
asymmetric networks. 

VIII. Appendix 

A. Proof of Lemma \3\ 

Define the largrangian of the sum rate maximization prob- 
lem for Tx i to be C{'Wi,iJ,) ~ R — /^maxlwf^Wj — 1). The 
neccessary condition of largrangian „ ^jj C{-Wi, /i) — gives: 



H i -Q^Rj — l^max^i- With elementary matrix 



calculus. 



9w 



H 



Ri 



P 



Efc=i|vfH,,w,pF + ai 



-E,w, (32) 



R, = - 



|vfH,,w,fP 



Ef^,|vfH,,w.pp + a| 



Erii|vfH,,WfcPP + a| 

-A,-W; 



(33) 



where A°f* is a function of all channel states information 



and beamformer feedback: 



Ea^i Sjk + CTj^ Y.k^j Sjk + a 



(34) 



where Sjk = |v^HjfcWfepP. Thus, the gradient is zero for 
any W; eigenvector of the matrix shown on the L.H.S. of 
( fTTb . Among all stable points, the global maximum of the cost 
function is reached by selecting the dominant eigenvector of 

B. Proof of Theorem \3][ convergence points of DBA 

To prove that interference alignment forms a convergence 
set of DBA, we will prove that if DBA achieves interference 
alignment, DBA will not deviate from the solution (stable 
point). 

Assumed interference alignment is reached and let 

lA „JA\ r- T aDL f^JA 



Let Qf^ = E^^.H 
Given receivers (v("^ , 



(w 



and (vf 



IA,,JA,HrTH 
"-tk 



and Q 



UL 



UL _ 



JA\ 



, v^' j, we compute new transport 
beamformers. In high SNR regime, Xji — > — oo and DBA gives 
w, = l/"""(Qf ^) (O. By (|29ll, Qf^ is low rank and thus 
Wi is in the null space of Q^^. In direct consequence, the 
conditions of interference alignment (l3Tl i are satisfied. Thus, 
(wi,...,WArJ elA"^. 



submitted to ffiEE TRANSACTIONS ON SIGNAL PROCESSING 



8 



Given transmitters (w("^, . . . , wj^), we compute new re- 
ceive beamformers. The receive beamformer is defined as = 
argmax^^^ — "^hqdl^. — Since Qf^ is low rank, the 

optimal Vj is in the null space of Qf ^. Hence, e IAF^. 

Since both and stays within lA^^ and lA^^ , 
interference ahgnment is a convergence point of DBA in high 
SNR. 



[23] G. N. He, M. Debbah, and S. Lasaulce, "K-player bayesian waterfiUing 
game for fading multiple access channels," in IEEE International Work- 
shop on Comutational Advances in Multi-Sensor Adaptive Processing, 
2009. 
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Fig. 5: Sum rate performance for asymmetric channel, with 
one link under strong noise, is illustrated. The strong noise, 
from out of cluster interference, is 20dB stronger than other 
links. DBA outperforms standard lA methods thanks to a 
proper balance between egoistic and altruistic beamforming 
algorithm. 



(b) Asymmetric uncontrolled interference and interference within cluster 
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(c) Asymmetric desired channel power 
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Fig. 4: Scenarios where asymmetry of channels are illustrated. 

Fig.l4alillustrates asymmetry of uncontrolled interference. Fig. Fig. 6: Sum rate performance for asymmetric channel, with 
l4blillusti-ates asymmetry of both uncontrolled interference and one link under strong interference within the cooperating 
interference within cluster. Fig |4c] illustrates the asymmetric cluster, is illustrated, 
strength of desired channel power 
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Fig. 7: Sum rate performance for asymmetric channel is 
illustrated. The direct channel gain of link 1 is 30dB weaker 
than other Unks. 



